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CSJ Abstract 

^ Quantum Portfolios of quantum algorithms encoded on qbits have re- 

cently been reported. In this paper a discussion of the continuous variables 

O version of quantum portfolios is presented. A risk neutral valuation model 
for options dependent on the measured values of the observables, analo- 
(~| gous to the traditional Black-Scholes valuation model, is obtained from 

f the underlying stochastic equations. The quantum algorithms are here 

I encoded on simple harmonic oscillator (SHO) states, and a Fokker-Planck 

, ^ I equation for the Glauber P-representation is obtained as a starting point 

for the analysis. A discussion of the observation of the polarization of a 
y—^ portfolio of qbits is also obtained and the resultant Fokker-Planck equa- 

^ tion is used to obtain the risk neutral valuation of the qbit polarization 

portfolio. 

00 

O 1 Introduction 

^' 

Recently, quantum portfolios of quantum algorithms encoded on two-level sys- 
terns (qbits) have been reported [IJ. It has been shown that maximizing the 
[ efRciency of quantum calculations can be accomplished via the formation of 

^ portfolios of the algorithms and by minimizing both the running time and its 

uncertainty. Furthermore, these portfolios have been shown to outperform single 
algorithms when applied to certain types of algorithms with variable running 
time such as the "Las Vegas" algorithms and NP-complete problems such as 
combinatorial search algorithms. In this paper a discussion of quantum portfo- 
lios of algorithms encoded on qbits and their continuous variables (CV) version 
(here, harmonic oscillators) is presented. Continuous variables quantum compu- 
tation performed with harmonic oscillators is discussed in [31 [31 . A quantum 
portfolio of algorithms is formed and a risk neutral model, analogous to the 
traditional Black-Scholes and Merton [3 (H [7] options valuation model, is ob- 
tained from the underlying stochastic equations. The quantum algorithms are 
here encoded on simple harmonic oscillator (SHO) states, and a Fokker-Planck 
equation for the Glauber P-representation is obtained as a starting point for 
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the analysis. A portfolio of qbits is also formed and the resultant Fokker-Planck 
equation of the qbit polarization is used to obtain the risk neutral valuation 
of the portfolio and measurement option. The results should prove useful in 
quantum computation and decoherence. 



2 Preliminaries 

A simplified model for quantum computation is proposed wherein the algorithms 
are encoded on simple harmonic oscillator basis states and are in the presence 
of a thermal bath, the temperature decoherence effects needing to be taken 
into account as is usually the case in real-world applications. The quantum 
computation utilizing harmonic oscillator basis sets representations has been 
discussed elsewhere [21 [3] . Representations of computation operators or effective 
operators can be made and the basis sets of these operators can be expressed as 
harmonic oscillator basis sets. For this case, we will also consider damping, this 
could be from other noise sources, external or system specific. The harmonic 
oscillator can be described by a master equation such as 
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Here, a and are destruction and creation operators as usual, and p is the 
density matrix operator. The mean thermal excitation due to the thermal bath 

is n = (e'=B^ — 1)~ , and F is the damping rate [8]. To obtain a Fokker-Planck 
equation for the SHO, we first suppose that p has a Glauber P-respresentation 

E 

^ (2) 

Substituting Eq.([2| into Eq.([T]) we obtain the Fokker-Planck equation with 
mixed derivatives 

(3) 

The Fokker-Planck equation is put into a regular form if we change variables 
using quadratures a = x -\- iy^ and Eq.([3]) becomes 

(4) 
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3 Risk Neutral Valuation 



The Fokker-Planck equation for the SHO impUes that there is an underlying 
stochastic differential equation(s) of the Ito form [8] 




fix 



dy= ^,fdt + ^ ^dW-,. 



(5) 



(6) 



Here, the Wiener processes dW (noise) are taken to be a Gaussian white 
noise, and are delta correlated. These stochastic processes can now be included 
in a portfolio analysis. To construct a portfolio, we first look at the case where 
the algorithm is encoded on two harmonic oscillators. It can be seen that a 
generalization to the multi-asset case proceeds straightforwardly. We write the 
Legendre transform for the N-asset portfolio 11 

A A 
n(*) = fixi...X,\,J/i...Tj,\.,t) [^A^J.; + ^,/J}l].. 

(7) 

and here, with N — 2 , the function /(xi, X2, 1/1^1/2, t) is the analogue to the (call) 
option in finance. In our case, it can represent the probability distribution for 
measurement of the observables of the quantum algorithm(s). These observables 
are Legendre transformed such that the state function, the portfolio 11 , evolves 
at the known or wanted rate r such that dll = rlldt. The N algorithms evolve 
independently, regardless of initial phase, and will have differing instantaneous 
values in their stochastic sample paths. They all have in common the same value 
for their moments (means and variances), as from the point of view of their 
statistics they have the same drift and diffusion coefficients and are equivalent 
on the average. This is different from the financial portfolio case where the 
portfolio of shares of a stock(s) assets have the same instantaneous value for 
all the shares (say, of a particular stock) and would resemble a portfolio of 
different traded stocks of one share per stock that are different in prices at any 
instantaneous period in time yet that equivalently have the same values for their 
1st and second moments of means and variances. 

Next we expand df in a Taylor series to first order in time, and to 2nd order 
in the variables, since dW^ = dt 

(8) 

Setting N — 2 for the two asset case, we take the differentials for dll. We have 
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(9) 

and we substitute the expressions for df and for dx'^s and dy'^s . Thus we are 
able to make the choice for the ^xi,yi niuhiphers directly [5j|6] as the terms. 

This Immediately eliminates the drift terms in the original process, in favor 
of the known rate r drift terms. This is a general feature in the sense that much 
more complex (nonlinear) process drifts can be eliminated in favor of the known 
rate r linear drifts [7]. Upon making the substitutions, we obtain the backwards 
Fokker-Planck equation that corresponds to the Black-Scholcs options pricing 
model. This PDE for the two-asset (and similarly for N>2 algorithms) portfolio 
is 



(10) 

The two-point version of this equation can be solved with the following choice 
for the function 



f[xi,Xi,yi,yi x[,x'2,y't,y'i,t') =e "fl(a:i,a;2,|!(i,s/2 x[,x'2,y[,y'j,t'). 



(11) 



With this choice, the PDE is transformed to a standard form backwards Fokker- 
Planck PDE. This backwards PDE is solved by g. Here the two-point function 
g defined in Eq.(ll) also solves the forward Fokker-Planck equation 

Thus the underlying stochastic processes are of the form 



(12) 



rfa- — rxdt + 



-dWi, 



(13) 



dy ■■ 



(14) 

The solution of the one-point function g is a Gaussian. We obtain the 
solution by transforming the stochastic differential equations to a drift-free 
dx => dzi,dy => dz2 form. The resultant one-point function f{zi,Z2,t) is 
then 



p 2D[e 



Anil 
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l]/(-2r)^ 



(15) 
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and similarly for the other variable z^. Here the diffusion constant is related 
to the damping rate as _D = The observables statistics 2-point function is 
formally solved by 

(16) 

and we can make the boundary value choices for the observables at the later 
time t' 

g{x\,x',,n].ji'.2,t") = S[x\ - Xi)5{x', - X.,)ii{y[ - Yi)Siy!, - ¥■,), 

(17) 

where these boundary values correspond to the terminal values of the evolution 
of the algorithms. These boundary values allow us to immediately integrate and 
obtain the solution for t' ~ T 

(18) 

The interpretation of this formal result for the delta function boundary values 
is that the value of the distribution of the observables of the quantum computa- 
tion algorithm is evaluated at the later time t' = T and then evolved backwards 
to the present time t similarly to the evaluation of the price of a financial in- 
strument such as the call or put option in portfolio theory. Other terminal 
boundary values can be made, such as step functions in direct analogy to insur- 
ance instruments. Closed form solutions are also possible, as the step function 
limits the integration and can be solved by several methods, and will not be 
discussed further here... note that due to risk management strategies in portfolio 
theory having a direct correspondence to observation strategies for portfolios 
of quantum computation algorithms, alternative boundary values are likely to 
become important for different applications. Also, as there are generalizations 
of the portfolio theory to 'beyond' the risk neutral valuation model, the simple 
delta function boundary values and formal solution can be considered as a sim- 
plest observation strategy . An example of the implementation of the strategy 
of observations of the quantum computations is the analog of the delta hedging 
strategy. The deltas ^ and ^ should be kept as close to zero as possible 
to ensure the evolution of the portfolio as dli = rlidt. This can be done by 
observing algorithms (selling) or adding algorithms (buying) to the portfolio at 
each computation time step. 
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4 Polarization of qbits and Fokker-Planck equa- 
tions 



In a paper that describes realistic models of interferometric continuous measure- 
ment of qbit readout systems relevant to quantum computing 9 , the author 
was able to describe the polarization of the qbit in terms of a Fokker-Planck 
PDE that accounted for the randomness inherent in the measurement process. 
The result due to the interferometric measurement Markov process was given 
by (a similar analysis can be performed for the qbits in the presence of thermal 
noise) 



here D{z) = (1 — z^)^/2 is the diffusion coefficient and V{z) — 2z{l — z^) is 
the drift velocity. The microscopic qbit polarization yields the time dependent 
polarization Z{t){—1 < z{t) < 1) since here 



(20) 



is the polarization from Eq. ( 20 1 . The other terms in Eq. ( 19 1 are the photon 



flux in photons per second, and 9 {order (10^ 6)) , the phase shift due to per- 
turbation of the qbit by a photon. The underlying Ito stochastic differential 
equation can be readily obtained and is 



(21) 



and can be seen to contain complicated drift terms as well as the diffusion coeffi- 
cients. This however will be the replacement by a linear drift term mentioned in 
the previous derivation. The portfolio analysis proceeds as before and we obtain 
for a two-qbit measurement distribution the following Black-Scholes PDE 



1=1 



(22) 



Defining the function / as / = e ^^'^ and substituting, we obtain the PDE 

2 



i—l 



(23) 
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This equation implies that the underlying stochastic process , say for dzi, is 
now of the form 



(24) 

The solution of the one-point function g is (for one asset) again a Gaussian 
form. The drift-free transformation xi — e~''*zi and X2 = e~'"'z2 allows us to 
obtain the solution for the one-point function g{xi, t),i = 1 



p SOKc-^" -ll/(-2rJJ(l-i,2j2 



V2^L»[fe - l]/(-2r))[l - ;?:r)- 

(25) 

The choice of the boundary terms is again the simplest possible delta functions 
corresponding to the instantaneous measurement. Performing the integration 
we obtain the formal solution 

f{2„z.,,t) = ft '•{'' *^g{z„z.,^tZ,..Z.,,t' = T). 

(26) 



5 Conclusion 

In this article, we have derived a risk neutral valuation model analogous to the 
Black-Scholes and Merton derivatives valuation model. This model is derived 
for portfolios of assets of quantum algorithms in continuous variables. The 
Fokker-Planck equations for the qbits are the starting points for the derivation 
of the formal solutions for the derivatives which in this interpretation are the dis- 
tributions for the measurement of the observables of the quantum computation 
algorithms of qbits encoded on harmonic oscillators and optics, risk neutral rate 
of evolution of the value of the portfolio dll — rlldt and for which the statistics 
of the computation are derived for future terminal computation time t' = T 
formal solution brought back to the current time t at the known or desired rate 
of return. This can be interpreted in several ways, including r the probability 
of successful computation per computation time step. The number of successful 
computations algorithms per unit time step, etc. and is an analogy to the rate 
of return on investment in finance. A connection to financial risk analysis and 
quantum computation is thus established, as the traditional trading strategies 
of risk minimization by portfolio optimization hedging schemes such as delta 
neutral, etc., have a direct correspondence to observation strategies. 
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